Prom Cooper pairs to Luttinger liquids with bosonic atoms in optical lattices 
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We propose a scheme to observe typical strongly correlated fermionic phenomena with bosonic 
atoms in optical lattices. For different values of the sign and strength of the scattering lengths 
it is possible to reach a "superconducting" regime, where the system exhibits atomic pairing, or 
a Luttinger liquid behavior. We identify the range of parameters where these phenomena appear, 
illustrate our predictions with numerical calculations, and show how to detect the presence of pairing. 
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After the achievement of Bose-Einstein condensation 
with dilute atomic gases ||], a great deal of interest in 
atomic physics has turned to the theoretical and experi- 
mental study of cold fermionic atoms. In particular, one 
of the most challanging goals nowadays is the observa- 
tion of the BCS transition Q , where fermionic atoms are 
expected to form Cooper pairs when they experience an 
attractive interaction. Although Fermi degeneracy has 
already been observed in several labs ||, the tempera- 
tures at which Cooper pairs are formed have not been 
reached (see, however, 

In this letter we propose to use bosonic atoms trapped 
in optical lattices to observe typical fermionic correlation 
phenomena, such the formation of Cooper pairs or the 
characteristic spin-density separation corresponding to a 
Luttinger liquid. Our scheme is based on the fact that 
strongly interacting bosons in a lattice can behave as 
weak (or even strong) interacting fermions. Our proposal 
is motivated by the recent experiment with bosons in 
optical lattices ||, where the strong interaction regime 
has been achieved, as it was theoretically predicted ||. 
In fact, this experiment illustrates that our proposal can 
be implemented with existing technology. 

The fact that strongly interacting bosons can behave 
as fermions is, of course, not a new idea. It is well known 
that the problem of one dimensional hard core bosons is 
exactly mapped into the one of free fermions |Q, since 
the infinite on site repulsion plays the role of an effective 
Pauli principle for bosons. This is the case, for instance, 
for the Mott phase of bosonic atoms in an optical lattice 
|^|, [|, or for the Laughlin atomic liquids in rapidly ro- 
tating traps ||. In an attempt to go beyond the usual 
free-fermionic behavior of hard core bosons, we show that 
the ability of tuning the sign and strength of atomic in- 
teractions can easily make bosons behave as interacting 
fermions. To be more specific, we will consider bosonic 
atoms confined in an optical lattice in such a way that 
tunneling can only occur in one spatial dimension. We 
will assume that some internal levels have been shifted up 
in energy using off resonant microwave fields, so that only 
two internal levels can be occupied. We will concentrate 
in the situation in which the interaction between atoms 
in the same internal level, U, is repulsive and stronger 
than the interaction between atoms in different ones, V 
H. The basic mechanism behind our scheme is simple. 



At sufficient low temperatures and tunneling amplitude, 
the strong repulsion U prevents two atoms with the same 
internal state to be in the same lattice site, so that atoms 
behave as interacting fermions. We will show that the in- 
teraction between the effective fermions is just the bare 
interaction V, and depending on its sign and strength a 
rich spectrum of possibilities can be observed. 

We consider a gas of N bosonic atoms in a one- 
dimensional optical lattice with M wells, and will call 
v = N/2M the filling factor. We will denote by a = f, I 
the two relevant internal levels. The Hamiltonian de- 
scribing this situation is H: 
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Here a\ a and bosonic operators that create (anni- 

hilate) an atom on the z-th lattice site with spin state <r, 
and riia = a ia ai a . The tunneling amplitude, t, as well as 
U and V can be easily written in terms of the scattering 
lengths and lattice parameters ||. 

Let us analyze first the simple limiting case U — » oo, 
where the atoms with the same spin are not allowed to be 
at the same lattice site. The ground state and low-energy 
excitations of Hamiltonian (|l|) lie within the subspace 
generated by states of the form: 
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where iVj + iVj_ = N, x = (xi,... ,xn t ), and y = 
(j/i, . . . , VNi) denote, respectively, the sites occupied by 
the particles with spin up and spin down, and Xi ^ Xj, 
Di Uji 3- This projected bosonic Hilbert space 

is isomorphic to the Hilbert space of fermions with spin 
1/2 in a one-dimensional lattice. The bosonic operators 
can be transformed into fermionic ones by the well known 
Jordan- Wigner transformation |llj : 
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where the operator cj CT creates an effective fermion in 
the j-th site with spin a. With the transformation (||) 
the bosonic Hamiltonian ( HI) p rojected to the subspace 
(0) is exactly transformed [H]] into a Hubbard chain of 
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fermions: 

H F = -t ]T c l c J° + Vj2 c h c ^ c li^l, (4) 

<i,j>a i 

where the parameters t and V coincide with those ap- 
pearing in Hamiltonian ([!]). 

Given the fact that the projected Hamiltonian and 
([|) are the same, they have exactly the same spec- 
trum. The bosonic eigenstates are related to the 
fermionic eigenstates by the following correspondence: 

= Ex, y M x -y) l x >y)F Ex, y /s( x ,y) l x >y) = 

|* B >. Here, |x,y) F = 4 lT . . . . . . c^JO), 

|x,y) is given by (§), and / F ( x ,y), /s( x ,y), the 
fermionic and bosonic amplitudes for the configura- 
tion (x, y), are related by a sign factor in the form 
/ fl (x,y) = sgn(x,y)/ P (x ) y) = J]^ sgn(x i -x i )sgn(y J - 

%)M x ,y) ; 

The physics of the fermionic Hubbard chain has been 
extensively studied (see, e.g., p2[). Depending on the 
sign of V and on the ratio V/t the system is known to 
exhibit different phenomena. In order to determine which 
of these phenomena can be observed with a real bosonic 
system we have to address two questions. First, since U 
cannot be infinite, we have to determine the conditions 
on U such that the above treatment remains valid. Sec- 
ond, we have to analyze the effect of the boson-fermion 
transformation (^|) on the predicted phenomena. 

Let us first estimate under which conditions our treat- 
ment will remain valid for finite U. Since we are inter- 
ested in the phenomena related to the ground state (and 
may be low-energy excitations) we have to compare the 
energy of these excitations with the one related to leav- 
ing the projected bosonic subspace. The first one can 
be determined from the known results for fermions [ fl2"| , 
whereas the second one will be of the order of U. Thus, 
we arrive at the condition U <C min(i, |V|). 

Now, let us discuss which fermionic phenomena can 
be observed with the bosonic system, i.e. the effects of 
transformation (^). These phenomena are characterized 
by the nature of the excitation spectrum and by some 
special behavior of the correlation functions. In partic- 
ular, for V > and v ^ 1/2 we have a Luttinger liquid 
p3| where charge and spin excitations are independent, 
gapless, and with phononic dispersion relations, leading 
to two different velocities. For V > and v — 1/2 there 
is a gap for charge excitations and the system is a Mott 
insulator. For V < the system is a superconductor. 
Fermions get paired, opening a gap for spin excitations. 
As the strength of the attractive interaction increases, 
the system evolves continuously from cooperative Cooper 
pairing (the BCS regime) to independent bound-state 
formation (the BEC limit) which is reflected in the 
correlation functions. We now argue that all these phe- 
nomena can indeed be observed with bosonic systems as 
well, i) Gaps for spin and charge excitations. The spec- 
trum for the bosonic system is the same that for Hubbard 
fermions. Moreover, spin and charge fermionic excita- 



tions are mapped onto spin and charge bosonic excita- 
tions, since all density (spin) operators c\ a Ci a {c\,c^) are 

mapped onto bosonic operators, amebic (ali a ii), by the 
transformation (|§|). ii) Correlation functions. Since the 
fermionic and bosonic amplitudes are related by a sign 
factor, it follows that any measurement in our bosonic 
system involving densities will give the same result as 
in the fermionic system, since the sign factor is then 
squared. Thus, both density-density correlation func- 
tions {n x +&n x ) , (where n x = n x -\ + n x i), and spin-spin 
correlation functions (S x +aSx) j (where S x = n x -\ — n x i), 
remain unchaged by (0). 

Note that not all observables are invariant under the 
transformation (Q). For example, the sign factor differen- 
tiating bosons and fermions will show up in the one-body 
correlation function (a< x+Aa a xa -) 1 which will be different 
for fermions and bosons. We can say that we have a sys- 
tem of bosons that behaves in many ways as a system of 
Hubbard fermions, though the bosonic nature of the real 
components of the system is not completely hidden and 
can be detected with some particular measurements. 

The above considerations imply that a variety of phe- 
nomena could be observed for the bosonic system. For 
example, for V > one may use these bosonic systems to 
observe spin-density separation in the same way as pro- 
posed for fermions [|l5| by just exciting locally either the 
spin or the density with a laser and looking at the prop- 
agation. Alternatively, one can analyze the excitation 
spectrum to see the two phononic branches. For V = 0, 
if we include an "impurity" atom at some position which 
strongly interacts with the rest of the atoms the system 
will exhibit the Kondo effect [[l6) . 

In view of the current interest in the achievement of the 
BCS regime with ultracold atoms, the rest of the letter 
will be devoted to the attractive regime V < 0. Although 
the problem can be exactly solved via Bethe Ansatz [18| , 
a variational formulation proves to be useful. We take 

N/2 

/u(x,y) = sgn(x,y) Y[ <fi{xi - yi), (5) 
i=i 

where ip is to be determined by minimizing the energy. 
The physical meaning of (|^) can be made more trans- 
parent if we write it in terms of fermionic operators as 

l*f> = E? r iP(* - v) 4i4d N/2 \°) which has been 

used in the literature for the fermionic case H], and 
even found reasonable agreement with the exact wave- 
function in ID fl7| . This state describes a condensate of 
pairs with spatial wave function <p, A = x — y being the 
relative distance of the pair. 

Let us first consider the regime of strong attraction 
|V| ^> t. Here, the sign factor in (||) is the same for all 
configurations (x, y). The bosonic paired state takes then 
the simple form |* B ) oc P[E£Li 4t4|] W/2 I )> where V 
denotes projection onto the subspace spanned by (^). In 
this state atoms are bound into on-site triplet pairs, and 
the system can be visualized as formed of hard-core com- 



3 



posite bosons. As it happens for Hubbard fermions with 
strong attraction the composite bosons move only 
via virtual ionization with an effective hopping ampli- 
tude tb — 2t 2 /\V\ <C t. When the attraction is so strong 
that we approach the limit \V\ = U the system looses 
its fermionic character. The state made by on-site pairs 
becomes unstable with respect to having more than one 
pair on the same site, since the cost in repulsion is ex- 
actly compensated by the gain in attraction. In fact, for 
\V\ > U it is favorable to have all the particles (N/2 with 
spin |, N/2 with spin J.) on the same site. The system 
behaves as a single "big boson" that moves very slowly 
along the lattice. As the ration \ V\/t decreases, the func- 
tion ip evolves continuously from a localized state of the 
Cooper pairs at A = to a delocalized one. 





FIG. 1: Projection of the exact ground state onto the one 
corresponding to the limit U — * oo as a function of |V|/t/ for 
N = 6, M = 8, and different values of t/U. 

To illustrate these behaviors we present now numerical 
exact results for a system of N = 6 bosons and V < 0. 
Even though this is a small number, this may closely 
represent some experimental situations with atomic sys- 
tems |5j where the effective number of atoms in the ID 
lattice may be relatively small. Given that Hamiltonian 
(|l|) is invariant under global spin rotations, we diagonal- 
ize it within subspaces of fixed total spin, S, and fixed 
z component of the spin, S z . In order to to check the 
"fermionization scheme" we show in Fig. [l] the projec- 
tion of the exact ground state of the system on the one 
calculated with the projected Hamiltonian, as a function 
of the strength of the attraction, \V\. In agreement with 
the predictions, when t <C U and as far as \V\ < U the 
overlap is more than 95% . However, as we approximate 
the limit |V| = U, or as t is increased, the system looses 
the fermionic character. 

For a regime of parameters in which the fermioniza- 
tion is valid, we can use the variational wave function (|^) 
to understand the nature of the bosonic pairs. Figure 
U shows the probability density of the pair, \ip\ 2 , corre- 
sponding to the variational state (S). We clearly see how 



FIG. 2: M 2 as a function of A/U for the variational wave 
function (pi) and for N = 6, t = 0.1U and different values of 
\V\/U. 



the nature of the pair exhibits a smooth crossover from 
being delocalized, for |V| <C t to get localize at A = 
for \V\ 3> t. The gap for spin excitations, given by the 
expression S = (-Bjy-i - 2E N + E N+1 )/2 Q, is plotted 
in Fig. H as a function of \V\. As the attraction increases 
the energy required to break a pair increases. The be- 
havior of the gap near |V| = U is due to the transition 
from a paired state to the "big boson" state. 
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FIG. 3: Gap for spin excitations, 8/U, as a function of |V|/f/ 
for iV = 6, and t = 0.1C7. The three regimes are indicated in 
the figure. 



The existence of pairs and the smooth evolution from 
Cooper pairs to localized pairs can be detected by an- 
alyzing the response of the system to a (Raman) laser 
that couples the internal states of the atoms, so that 

H B = H + BYZ ■i a \] a il a li a it> where the parame- 
ter B can be varied by tuning the intensity of the laser. 
The effect of this laser can be understood as follows. Let 
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FIG. 4: <Sf) as a function of B/U, for N = 6, t = 0.1(7, and 
different values of V/U. The inset shows the second deriva- 
tive, 77, of (S%) with respect to B (in units of U) at zero 
magnetic field, as a function of |V|/[7. 



us assume that the system is in a paired state with pair 
wave function >p(x)(\ Tl) + | IT))- The laser will try to 
rotate the spin state into the state | Tt) + I ID- Since 
this rotation must break the pair it follows that the re- 
sponse will be very low until a value B ~ 5 is reached. 
Thus, 5 can be detected by observing, for instance, the 
fluctuations induced in the the z component of the to- 



tal spin, (5*2), as a function of B, as it is illustrated in 
Fig. |J. As predicted, the system response, quantified by 
d 2 {Sz)/dB 2 \ B=Q , is extremely large for V = 0, when we 
have no pairs and the system holds no resistance to be 
rotated, and decreases monotonically for increasing at- 
traction, revealing the existence of pairs. Note that for 
\V\ t, rotation of the on-site pair will imply exciting 
the system out of the subspace spanned by (g) , since two 
atoms with the same spin will be left on the same site, 
and therefore excitation will occur for B > U. 

In summary, we have shown that with the experiments 
that are being carried out with bosons in optical lattices 
one can observe a great variety of fermionic correlation 
phenomena, including Cooper pair formation and spin- 
density separation. These atomic systems may be simpler 
to manipulate than fermions themselves and can provide 
a very clean laboratory to study several interesting phe- 
nomena which have eluded observation so far. Apart 
from that, the system proposed here may bring up novel 
phenomena in 2 and 3D optical lattices. Even though the 
Jordan- Wigner transformation is not appropriate, some 
of the phenomena predicted for Fermions (like d-wave 
pairing or an analogue) may appear in the case of bosons 
in the regime U > \V\, and therefore these cases are 
worth exploring both theoretically and experimentally. 

Discussions with I. Bloch, J. von Delft, C. Tejedor, and 
P. Zoller are gratefully acknowledged. B. P. is indebted 
to G. Gomez-Santos for illuminating discussions. 
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